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Abstract 

We consider gravitino couplings in theories with broken gauge symmetries. In particular, we 
compute the single gravitino production cross section in W~^W~ fusion processes. Despite recent 
claims to the contrary, we show that this process is always subdominant to gluon fusion processes 
in the high energy limit. The full calculation is performed numerically; however, we give analytic 
expressions for the cross section in the supersymmetric and electroweak limits. We also confirm 
these results with the use of the effective theory of goldstino interactions. 



1 Introduction 



One of the reasons that supersymmetric theories are the prime focus for physics beyond 
the standard model, are their inherent abihty to be tested. In addition to its many more 
theoretical benefits such as stabilization of the electroweak symmetry breaking scale [1] and 
unification of gauge couplings at high energy [2] , low energy supersymmetry [3] often predicts 
a particle spectrum readily observable at colliders such as the LHC Indeed, in models 
where unification conditions are placed at some high energy scale, such as in the constrained 
minimal supersymmetric standard model (CMSSM), regions of parameter space consistent 
with known phenomenological constraints at the 95 % CL are well within the expected reach 
of the LHC [5]. If i?-parity is conserved, supersymmetry (SUSY) also predicts that the 
lightest supersymmetric particle (LSP) is stable, thereby making it an excellent candidate 
for dark matter and if it is the neutralino [6], it is also potentially observable in direct 
detection experiments [7]. In this case, it is usually assumed that the gravitino is heavier 
than the neutralino, and even then, additional assumptions must be made so that its decays 
in the early universe do not upset the results of big bang nucleosynthesis [8l|9]. 

It is also quite possible that the gravitino is the LSP [6l[T0l[Il]. In the CMSSM, this will 
occur whenever the gravitino mass, ^3/2 is less than the lightest standard model superpartner 
mass pAj making it subject to big bang nucleosynthesis constraints on the decays of the next 
to lightest supersymmetric particle (NLSP) [T^P^] . Indeed, a gravitino LSP is quite common 
in models based on minimal supergravity [H]. Typically, one would expect gravitino masses 
of order the weak scale, making direct detection of dark matter very unlikely. There are 
nevertheless proposals for detecting the long lived decays of a stau at the LHC p^llG] . 

The possibility of a light gravitino precedes the MSSM pTp8] and in models of gauge med- 
itated supersymmetry breaking [19], the gravitino may be significantly lighter with masses as 
low as 10"^ - 10"^ eV. While cosmological constraints on these models may be derived [20], 
there remains a broad mass range for super-light gravitinos. In no-scale supergravity mod- 
els [21] , the gravitino mass is decoupled from the rest of the supersymmetric sparticle spec- 
trum and may be set to the Planck scale [22], or to the keV scale and below [23]. 

The detection of very light gravitinos at colliders is in principle possible through the decay 
of the NLSP [211[25] or through direct production at e~^e~ [21] or hadron [25l[2I| colliders. 
This is possible, because, as we will see, the gravitino couplings are inversely proportional 
to its mass, making very light gravitinos more readily accessible. The current lower bound 
on the mass of a super-light gravitino comes from LEP and is [28] 




and this limit will surely be improved at the LHC. The dominant processes affected by a 
light gravitino are expected to proceed through the pair production of gluinos 




m3/2 > 1.35 X IQ-^eV, 



(1) 



PP 99 



(2) 



or through associated gravitino production with either squarks or gluinos 



pp^gG, qG 



(3) 
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It is known that gravitino production processes suffer a breakdown of unitarity at high 
energies due to the non-renormahzabihty of the super-gravity Lagrangian [29]. However, 
given the mass bound ([T]), unitarity is preserved through the TeV scale. Recently, it was 
claimed [30j that the breakdown of unitarity is significantly more severe in theories with a 
broken gauge symmetry such as the Standard Model. Indeed, it was claimed that in the high 
energy limit, the cross section for gravitino production remains non-zero even in the limit 
of exact supersymmetry. If true, this would imply that associated production of gravitinos 
through W boson fusion would come to dominate at high energy when compared to gluon 
fusion (where the gauge symmetry is unbroken). Here, we will calculate the weak boson 
fusion process leading to gravitino production and show that, contrary to the claims of [30] , 
the cross section is well-behaved at high energy. 

For very light gravitino masses, couplings of the gravitino to matter are dominated by 
the goldstino and gluon fusion process will therefore be proportional to m^/rn^i^- Single 

gravitino production through gluon fusion, pp — )■ gG was recently reconsidered in [31j, where 
they found 

assuming mass spectra corresponding to SPS benchmark points 7 and 8 [32] for which the 
gluino mass is 920 and 810 GeV respectively. 

In contrast, neutralino pair production through W^W~ fusion, pp — )■ x^'x"? was consid- 
ered in Ref. [SB]. Cross sections for the same SPS benchmark points 7 and 8 were found to 
be as high as ~ 0.1 fb (for the production of y^xV)- would naively estimate that single 
gravitino production would scale as 

where msusy is the typical mass of supersymmetric particles (close to the electroweak scale) 
and Mp is the reduced Planck mass Mp = ^ 2.4 x lO-*^*^ GeV , where Gat is the 

Newton gravitational constant. 

The ratio between the gluon and W fusion productions are therefore estimated to be 

^ ~ 10-^ (6) 
If the claim in [30] was right, one would get an additional factor of s/ml^^y-. 

Iperrantelli ~ 10 ^ ( ^ ) (7) 

where s is the square of the center of mass energy in the collision between W^W^ . This ratio 
can in principle be of order one, or bigger. In Fig. [H we show the qualitative behavior of 
the cross section for single gravitino production in the symmetric case of gluon fusion (solid 
curve labelled gg), and in the broken case of W boson fusion (dashed curve labelled WW). 
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If eq. ([7]) holds, W fusion process would come to dominate over gluonic ones (as shown by 
the dotted curve). Note that the cross sections shown in the figure do not include necessary 
form factors and so do not represent p p cross sections. 
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Figure 1: Single gravitino production cross section as a function of center of mass energy. 
The gluon fusion process leading to gravitino plus gluino is shown by the solid curve labelled 
gg. The cross section for W fusion to gravitino plus neutralino predicted in [30] is shown by 
the dotted curve. Our calculation for the same process is shown by the dashed curve labelled 
WW. Choices for the supersymmetric parameters used are discussed in section [31 

Furthermore, we point out that the claim of [30] is also in contradiction with the equiva- 
lence theorem, which states that the gravitino can be effectively replaced by the goldstino pT] 
at energies much greater than its mass. In calculations based on the equivalence theorem, 
one uses the on-shell conservation of the supercurrent, to which the goldstino is coupled to. 
If the result of [20] was correct, it would imply that a broken gauge theory contains some 
loopholes that invalidate the theorem, and that would not allow one to use the equations of 
motion in determining the couplings of the goldstino. This prompted us to revisit this issue, 
and to provide an explicit proof of the theorem. We do so by generalizing the calculation 
of [3l|, where the equivalence is shown at the level of 5— matrix elements. The calculation 
of [31] is performed for an unbroken U(l) theory, and the effective Lagrangian included a 
term which was previously unnoticed. The result of [31] is confirmed by [35] and the effective 
Lagrangian is related to the soft SUSY breaking terms in the MSSM through an explicit use 
of the equations of motion. Here, we extend these derivations to include the case of broken 
gauge symmetries. We make this relation explicit in equations (fTTjl and ( 125|1 below, and 
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WW 



prove it in Appendix [Bl It is manifest from the proof that the theorem apphes irrespectively 
to whether the gauge symmetry is or is not broken. 

The plan of the paper is as follows. In the next section, we write out the interaction 
Lagrangian for the gravitino coupled to the MSSM with broken electroweak symmetry. In 
anticipation of taking the high energy limit, where we can replace the gravitino couplings with 
the goldstino, we write out the explicit couplings of the goldstino to the terms originating 
from the soft supersymmetry breaking Lagrangian. In section [3l we outline our calculation 
of the W~^W~ — )■ Gx^ cross section. While the analytic expression for the cross section is 
too long to write out, we do give analytical results in a couple of interesting limits. First, we 
show that in the supersymmetric limit {7713/2 <^ ^«susy ^ ttiw), |A^P oc mg^gyS/mg^g^p at 
high energy which would lead to a cross section for gravitino production of the form given in 
eq. ([5]). We also consider the limit 7713/2 ^ niw ^ "^susy and write out the analytical cross 
section at high energy, which takes a similar form. We also comment on the detectability of 
gravitinos through this process in comparison with that of gluon fusion. Concluding remarks 
are given in section m We also show explicitly our derivation of effective gravitino Lagrangian 
in the appendix. 

2 Interaction Lagrangian for the gravitino with broken 
electroweak symmetry 

The interactions vertices between a single gravitino and the MSSM fields are obtained from 
the interaction Lagrangian 

^ _ 



V2M, 



8M, 

where, following the notation of [36], ^/^^ denotes the gravitino field, (j) and xl the scalar 
and fermion components of the chiral MSSM superfields, Fpcr is the field strength of a gauge 
boson field, and A is the corresponding gaugino. The indices i and (a) label the chiral and 
gauge multiplets, respectively (notice that we are implicitly summing over all the MSSM 
chiral and gauge multiplets). The covariant derivative of a scalar field is 

^ d^cj)^ + I 4")^^ (T(")'^ (f)' (9) 

In the first line of (|8]), S^ssm denotes the contribution from the MSSM fields to the 
supercurrent and contains only terms from the supersymmetric Lagrangian. Specifically, 
under a supersymmetry transformation, any MSSM field (of any spin) $j transforms as 
$j — )■ $j + while the supersymmetric part of the MSSM Lagrangian transforms as 
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Ausy — ^ -Csusy + d^K^. Theii, the supercurrent is 

dC. 



Ql^ = SUSy 

^MSSM - 5(5^$^)^^^ ^ ^^Uj 

The exphcit expression for the supercurrent can be found for example in [35] . 

We want to single out the gravitino interactions that arise due to the breaking of the 
electroweak symmetry. Following [37], we denote the two Higgs doublets as 

■ "'^[m) ' '''' 

and we denote their vacuum expectation values (vevs) as {Hi) = vi , = V2 , {Hf) = 

{Hi) = 0. We denote the corresponding higgsino fields as 

HIl = PlH^ , HI^ = PlH~ , HIl = PlH+ , HIl = PlH2 , (12) 

where H^ = (j^^ 

Then the interaction term we are interested in is 



Ant 3 



vi W^^ PlH- + V2 W^- PlH^ 



f^= — (vi PlHi - V2 PlH, 

yl cosy-.. ^ 



+ h.c. (13) 



(where we have used ipij.'~i^ = 0, and where the gauge fields are defined in the standard way, 
see for instance [3j). We rewrite these interactions in terms of the chargino {xj ,J = 1, 2) 
and neutralino {x^ ,i = 1, . . . ,4) mass eigenstates, using the rotation formulae [37] 
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PlH^ = V*2PlXj , PlH- = U;2PlX'j 
PlH, = N*,Plx- , PlH2 = N*,Plx', . (14) 

We also rewrite the two Higgs vevs in terms of the (tree level) M\y and Mz, using the 
notation [37]: = {vf + f^) /2 and tan/3 = V2/V1. We end up with 

Ant 3 -^tP^PL[V2Mw{cos/3U;2W^+x- + sm(3V;2W^-Xj) 

+Mz (cos /3 iV*3 - sin /3 N*^) xf\ + h.c. (15) 

The remaining interactions between the gravitino and MSSM fields coming from (|8]), 
can be found in [36j using MSSM gauge eigenstates in the absence of electroweak symmetry 
breaking. In AppendixjXl we rewrite the gravitino-MSSM interactions in terms of the MSSM 
mass eigenstates, including the effects of electroweak symmetry breaking in the rotation 
matrices (between gauge and mass eigenstates). 

The couplings of the gravitino at energies much greater than its mass can be more easily 
written in terms of an effective interaction between matter and the goldstino field [17]. The 
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situation is analogous to what happens for spontaneously broken gauge theories, for which 
the couplings of the longitudinal polarization of massive gauge bosons are determined at 
high energies by those of the goldstone bosons that are eliminated in the unitary gauge (this 
is known as the equivalence theorem). Analogously, the gravitino is coupled at the quadratic 
level with the goldstino field. In the super-Higgs mechanism, the goldstino is absorbed 
into a redefined gravitino field (or, equivalently, it is set to zero in the unitary gauge). 
At energies greater than the gravitino mass, the longitudinal gravitino component is more 
strongly coupled to matter than the transverse modes, and the couplings are determined by 
those of the absorbed goldstino field (for a recent general study of the phenomenology of a 
strongly coupled glodstino, see j38j). 

We can see this from the polarization tensor [39] 



(r) 




^ ' 7m + I (#-"^3/2) (7^^+ 



3 V ^3/2/ ^ V ^3/2 



(16) 



where p and m3/2 are the gravitino momentum and mass, respectively. In the last expression, 
we have written the leading term in the polarization tensor in a 1/777.3/2 expansion. This term, 
which comes from the longitudinal gravitino polarization, dominates at energies greater 
than the gravitino mass. Since X]r=±i X^"^^ X^"^^ = f + ^3/2 ~ f-, we can effectively replace 

\l\'mY '-^P irrelevant phase) in this high energy regime. Therefore (after an 

integration by parts) 

7 — 7 

^'^"^^ ^ „^ '5'msSM + h. C. — £int, eff = ^= — — X 'S'^SSM + h. c. (17) 

v2Mp v3m3/2Mp 

We can actually simplify this expression further, and obtain an effective interaction La- 
grangian in non-derivative form. To see this, consider the infinitesimal variation of the MSSM 
Lagrangian £mssm = -^susy + >Csoft under an arbitrary infinitesimal variation of the MSSM 
fields. Since only MSSM fields or their first derivatives enter in the Lagrangian, one has 

5£mssm = 5$. + 6d,^.. (18) 



One can then immediately rewrite this expression as 



- O, 



9$, ^ d {d^^i) 



Let us now specify the infinitesimal variations to be the variations of the MSSM fields 
under a supersymmetry transformation. Using eq. (fTOjl . and the fact that £soft does not 
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contain first derivatives of fields, we have 



dC 



MSSM 



5$i 



d.. 



dC. 



susy 



(20) 



where we recall that (9^ K^^ is the variation of £susy under an infinitesimal supersymmetry 
transformation. Therefore 



5C 



MSSM 



5C 



soft 



5C 



MSSM 



Inserting this equation in eq. fl20|) . and the resulting expression in eq. f|T9l) . we obtain 



5C 



MSSM 



dC 



MSSM 



d, 



dC 



MSSM 



or 



^MSSM 



dCsoit 
5$, 



-MSSM 



f^-^MSSM o <9>Cmssm 
- o, 



(21) 



(22) 



(23) 



Inserting this expression in eq. f[T7jl . we rewrite the interaction Lagrangian between the 
MSSM and the light gravitino as 



-"int, efF 



MSSM 
9$, 



9, 



MSSM 



+ h. c. 



(24) 



As we prove in Appendix [HI the part in square parenthesis does not contribute to the 
amplitudes of physical processes having one light gravitino in the initial or final state (in 
short, one can take the on shell expression for 5'^ssm; since the term in square parenthesis 
vanishes on shell; notice that the procedure just outlined provides the on-shell expression 
of 'S'^ssM without the need to explicitly work out the equations of motion of the fields 
entering in the supercurrent). Namely: 



int, cfT 



X — + h. c. 



(25) 



This is the effective theory for the MSSM-light gravitino interaction in non-derivative form. 
To get an explicit expression, we recall the MSSM superpotential and soft supersymmetry 
breaking Lagrangian: 

W = huH2Qu' + hdHiQd' + KHiLe" + IJH2H1 (26) 



soft 



- ( -M,A"A" + h.c. 



mirct^' (27) 
- {AuKH2Qu'' + AdhdHiQd" + AeKHiLe" + B^iH2Hi + h.c.) 
where generation indices on the matter fields have been suppressed. From this, we find 



int, efT 



4 \/QMp ms/2 



a/3 Mpm3/2 



y/6 Mp m^/r 



(r7;;;.0^-)x7'A(")" (28) 



where {(pi, Xl) are the MSSM chiral multiplets, mjj are the (low energy) scalar mass^ terms. 
In the second term, Aj = A^, A^, A^, and _B for j = 1 — 4, and Wj refers to the respective 
term in the superpotential. The indices j and i are both summed and the latter runs over 
the chiral fields in each Wj. 



3 Computation of W — > gravitino + neutralino 

We now compute the cross section for the scattering of unpolarized W pairs 

W+ (k) + W- ik') ^Gip) + x° (g) , (29) 

which was studied in ^U\. There are five diagrams contributing to this process, which we 
show in Figure El The contribution from each diagram to the amplitude of the process can 
be found in Appendix [Cl 




Figure 2: Diagrams contributing to W'^ W — )■ gravitino + neutralino. 

We denote by A^j the matrix element for each neutralino mass eigenstate produced 
in the scattering. We computed the unpolarized squared matrix element 

spins, polarizations 

with the aid of the Mathematica package FeynCalc [31] (note that the index i is not summed 
over). The resulting exact expressions are too long to be reported here. 

The matrix elements are dimensionless, and since all gravitino vertices are proportional 
to 1/Mp, all of the squares of the matrix elements are proportional to 1/Mp. The sum over 
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the gravitino polarizations produces a term proportional to 1/7713^2, cf. eq. (fT6|) . As a con- 
sequence, the strength of the gravitino interactions increase with decreasing 7/13/2, and one 
can use accelerator constraints to set a lower limit on the gravitino mass as discussed earlier. 
In Ref. |30], it was claimed that the term proportional to l/mg^g results in a contribution 

l-^l? = M^^'^ unpolarized squared matrix elements, with a coefficient X that does 

not vanish in the limit of exact supersymmetry. More precisely, X was found to be propor- 
tional to the sum of a product of two Mandelstam variables. This is quite different from the 
case of gravitino production from gluon fusion, where the analog of X is proportional to the 
square of the gluino mass and hence vanishes in the limit of exact SUSY. The discrepancy 
was attributed in [30] to the breaking of the electroweak gauge symmetry. 

Although the complete result for the amplitude is too long to be reported here, we show 
in Fig. [T]the resulting cross section for a specific choice of parameters (see below for details). 
Moreover, we analytically study, and present, the result in two relevant limits. Several 
combinations of the physical masses and rotation parameters 

M^o, M^o, /3, a, M^^., Uij, f/2,, Vij, V2,, M^o, Na, N,2, Ni^, Ni^ (31) 



appear in our (very lengthy) exact results for |A^|f. a is the mixing angle in the scalar 
Higgs mass matrix, and the other quantities have been defined above. As a consequence, the 
behavior of the exact expressions in the limit of exact supersymmetry is not manifest. The 
supersymmetry breaking parameters in the MSSM that are relevant for this computation are 
the soft masses Mi, M2, rn?Hu^ "^l^d' ^ introduced in eq. (125]) . To study the supersymmetric 
limit of \M.\1, we assume that the soft SUSY masses are of the same order of magnitude, 
which we denote as msusy: 

Ml ~ M2 ~ |mH„| ~ Im^J ~ |5| = O (msusy) , m3/2 < msusy < Myy (32) 

The limit ^3/2 ^ ^susy is the phenomenological relevant limit for setting accelerator bounds; 
the limit mgusy ^ Mvi/, although not physically realized in Nature, is the appropriate as- 
sumption to analytically study the claim of [311]. 

In the limit 7/13/2 ^ ^susy ^ our averaged amplitudes can be formally written in 
the form 

where the coefficients of the expansions, (some of which may vanish) are independent 
of r/igusy and m3/2. Note that the coefficients will be different for the different outgoing 
neutralinos. 

The only other relevant input parameter besides the soft masses (1321) is the /x parameter 
of the Higgs potential. The minimization of the Higgs potential leads to two equations which 



^We stress that all the terms of the sum ([T6| are included in our exact expressions from which the 
expansion (|33| is performed. 
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allow one to solve for the two expectations values or equivalently, Mz = {g^ + g''^){v\ + v^) j 2 
and tan/3. Instead, it is common to specify Mz and tan/3, in which case it is possible to 
solve for the Higgs mass mixing parameter, /i, and -B, 



m 



m]j tan^ (3 + |M|(1 - tan^ /3) + A 



(1) 



(2) 



tan2 + 
i(m2 + + 2/i') sin 2/3 + As, 



(34) 



("1 2) 

where A^ and A/(' are one-loop corrections to and but will be ignored in our 
analytic expansions as we are restricting our calculation of the cross section to tree level. In 
the supersymmetric limit yU, i? — )■ and tan /3 — )■ 1 (see e.g. [37]). From these expressions, 
we can express the parameters m as expansion series in m,^,y/Mw. H Finally, we insert 
these expressions into |A^|^, and we find that the result 



Xoo = Xoi = Xio = in eq. m 



(35) 



is indeed recovered for all i = 1,2,3,4 (this explicitly shows that the numerator of fl33|) 
vanishes for exact supersymmetry, in contrast to what claimed in [30]). The remaining 
terms are in general nonvanishing. In the limit m3/2 <^ ^susy M]^, the unpolarized 
squared matrix elements are dominated by the term proportional to X20 in eq. fl5^ . Since 
X20 has mass dimension 2, and since the external momenta can be expressed in terms of 
Mandelstam variables, we formally have X20 = do x momenta^ + (^2 x M^, where (io,2 
are combinations of dimensionless quantities (in practice, only numerical factors and the 
gauge group charges). The full expressions for X20 are still too lengthy to be written here. 
However, the terms proportional to do, dominate in the high energy limit, and in this case 
we can write out analytic expressions for the unpolarized squared matrix elements. For 



\M\l 



\M\l 



\M\ic^\M\ 



eh M| 



27 sin^ Ml ml^^ 



s 



Mi(Mi + 2M2) (l + 2^)' + M| (l7 + 36f + 36f 



108 M2m2/2 



s 



216M|,m2/2 



X 



Mi(tan2^^Mi-2M2) (1 + 2- + 

s J sin f/^( 



4 + I 17 + 36 - + 36 — 



(36) 



cos 6',, 



In this limit, the masses of x\ and x\ go to (as m^usj ~^ 0). These states are a symmetric 
combination of the Higgsinos and the photino respectively. The masses of Xs and xa both 



^The resulting expressions for ([5T|) are lengthy, and we do not report them here. 
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approach Mz, and these are mixtures of the zino and an anti-symmetric and symmetric 
combination of the Higgsinos. 

We also computed the scattering W~ — )■ goldstino + neutrahno using the effective 
theory ( l28l) . and we precisely recovered the expressions (1361) in the high energy limit. We 
note that the first and last term in (|28ll also give quadratic goldstino-neutralino interactions, 
proportional to the two Higgs vevs. Such terms are included in the computation as mass 
insertions. 

It is also interesting to study the exact results in the limit of Mw -C ^susy, since this is 
the more phenomenologically relevant one. Repeating the same exercise discussed above, we 
find in the hmit 7713/2 <C Mw ^ m^nsy ^ a/N' 



\M\l 



\M\i 



108 cos2 Ml m\i^ 

108 sin^ Ml m\i^ 
e^Ml 



\M\ 



+ s cos^ (2/3) 



9 (s + 2ty 



+ 8 s + s cos^ (2/3) 



^ 27 sin^ Ml m\i^ 



(37) 



As one would expect, we again find that the numerators vanish for mgusy — ?■ in this limit. 
Now, the four neutralinos have masses which approach Mi, M2, [i and /i respectively and 
are effectively the bino, wino, and antisymmetric and symmetric Higgsinos. As we did for 
eqs. ( 136|) . we also reproduced the results (1371) using the effective goldstino Lagrangian ( 128|) . 

From both ( l36l) and ( 1371) we see that the square amplitude grows linearly with the Man- 
delstam variables. The resulting cross section is therefore constant at s^/^ much greater than 
the masses of the particles involved in the scattering. This is in contrast with the cr cx s 
dependence claimed in [30] . 

For illustrative purposes, we show the cross section for a specific choice of parameters. 
We choose to work in the context of no-scale supergravity [21] characterized by the Kahler 
potential 



G 
K 



K 

Ml 



+ F {(t^) + Ft (0*) 



-3 Ml In 



z + z* 
Mp 



3 Ml 



(38) 



where for simplicity we consider only one hidden sector complex field, z. The scalar potential 
takes a globally supersymmetric form 



(39) 



plus Z)-terms. It is important to note here the absence of all of the soft supersymmetry 
breaking masses. That is, at the scale at which supergravity is broken (which we assume to 
be greater than the grand unified scale), t/Jq = Aq = Bq = 0. These terms will be generated 
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radiatively from the non-zero gaugino mass, which at the supersymmetry breaking scale is 
given by 

1 



m\ = - 
2 



ef?/2^(lni?e/i): 



(40) 



where h{z) is the gauge kinetic function assumed to be diagonal in its gauge indices. For 
the no-scale Kahler potential, one then finds that 

For a suitable choice of h p3], the gravitino mass can be made much smaller than the gaugino 
mass. 

Phenomenological models based on no-scale supergravity have been recently constructed 
[l2], and we use two examples of low energy spectra based on that work. In the first example, 
we choose a supersymmetry breaking scale of Mp, and a universal gaugino mass mi/2 = 600 
GeV. Recall that mo = Aq = Bq = 0. The low energy spectra also depend on two couplings 
in the GUT scale superpotential corresponding to the term cubic in the Higgs adjoint (A') 
and a mixing term between the adjoint and the Higgs 5-plets (A). In this example, we take 
A = —0.06 and A' = 1. Because we are specifying Bq at the input scale, we are not free to 
choose tan /3. In this example, it is calculated to be tan /3 = 47.8. When run to the weak 
scale, this model has gaugino masses of Mi = 275 GeV, M2 = 534 GeV. The soft Higgs 
masses are Tn'j^^ = —1024^ GeV = —615^ GeV. When loop corrections are included 

in calculating the low energy spectrum, we find fi = 840 GeV, and neutralino masses of 
283,550,913 and 919 GeV. The gluino mass is 1510 GeV. The scalar Higgs masses are 119 
and 734 GeV. We have fixed the gravitino mass to 1713/2 = 10~^eV. 

We show in Figure [3] the cross sections for the production of a gravitino and each of 
the neutralino eigenstates. These cross section are evaluated numerically from the exact 
square amplitudes. We see that they indeed approach a constant value at high a/s. For this 
choice of parameters, the processes producing the first two neutralinos have a resonance at 
^/s = Mfjo , corresponding to the heavy Higgs exchange process. The resonance is narrow as 
compared to the range of ^/s shown here, and it is barely visible in the result for Xi shown 
in the Figure. 

Figure [3] shows the cross sections until a/s becomes too large, and our results are affected 
by numerical inaccuracies. However, one can verify that the analytic approximations written 
above are in excellent agreement with the exact expressions in the high energy limit. In Fig. 
m we compare the exact results with the approximations ( l37j) . for the illustrative case shown 
in the left panel of Fig. [3] where the soft supersymmetric masses are all sufficiently greater 
(in magnitude) than rriw In Figure HI we present the comparison for the process producing 
the second neutralino eigenstate (which is the dominant one for this choice of parameters). 
An equally excellent agreement is also found for the processes producing the other three 
neutralinos. 

To see the dependence on our particular choice of low energy spectrum, we show in the 
right panel of Fig. |3] a second example. In this case, we choose a supersymmetry breaking 
scale of 3 X 10^^ GeV, and a universal gaugino mass mi/2 = 300 GeV, with A = 0.2 and 
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Figure 3: Cross section for W~ — )■ gravitino + neutralino for the two sets of parameter 
choices specified in the main text. The gravitino mass has been fixed to m3/2 = 10~^eV. 
We recall that the cross section scales 
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Figure 4: Comparison between the exact cross section shown in the previous Figure (for the 
second neutralino eigenstate), and the one obtained from the approximated square amplitude 

m 



X' = 1. We find tan/3 = 11.2 and gaugino masses of Mi = 123 GeV, M2 = 239 GeV The 
soft Higgs masses are m% = —477^ GeV m% = 189^ GeV, /x = 413 GeV, and neutralino 
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masses of 123, 231, 426 and 442 GeV. The gluino mass is 745 GeV. The scalar Higgs masses 
are 113 and 464 GeV. Again, we have fixed the gravitino mass to ^3/2 = 10~^eV. 

In our second example, the sparticle spectrum is somewhat lighter (by roughly a factor 
of 2). As one can see, the qualitative behavior of the cross sections is similar to that found 
in the left panel. The heavy Higgs resonance however, is now much more prominent. 

4 Summary 

It is quite possible that low energy realizations of supersymmetry yields a spectrum with a 
gravitino LSP. In models with gauge mediated supersymmetry breaking, as well as in no-scale 
supergravity models, the gravitino may in fact be very light compared with the rest of the 
superpartner spectrum. As we have discussed above, and shown rigorously in Appendix B, 
the high energy interactions of a light gravitino are dominated by its longitudinal component, 
or goldstino. As a result, the couplings of gravitinos to matter are proportional to 1/ Mpm^i^ 
making them readily accessible in accelerator searches. 

One might expect that gravitino production cross sections mp p collisions be dominated 
by quark and gluon fusion process, however, a recent calculation claimed that due to effects 
associated with electroweak symmetry breaking, W boson fusion process would eventually 
come to dominate the overall gravitino production cross section at high energy. Here, we 
have shown this claim to be untrue. We have calculated the gravitino production cross 
section in both the high energy and supersymmetric limits and found no enhancement due 
to electroweak symmetry breaking effects. 

Although the full analytic expression for the cross section for gravitino production through 
W fusion is too lengthy to write out, we were able to express the cross section in two limiting 
cases both at high energy: the supersymmetric limit - where one sees explicitly the fact that 
the matrix elements are proportional to the gaugino masses (just as they are for the case of 
an unbroken gauge symmetry) and in the (more physical) electroweak limit where again the 
matrix elements are proportional to gaugino masses. 

We have also worked out in detail the applicability and use the equivalence theorem. 
In the high energy limit, the interactions of the gravitino can be replaced with derivative 
interactions of the goldstino. After an integration by parts, the goldstino is coupled to the 
divergence of the supercurrent. In Appendix [Bl we prove the equivalence theorem and show 
that its validity does not require an unbroken gauge symmetry. As a consequence, we are 
able to write down a relatively simple form for the effective interaction Lagrangian, and 
verify that the resulting cross section agrees with the original result. 
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A Single gravitino-MSSM vertices with broken elec- 
troweak symmetry 

We write here all the interactions of a single on-shell gravitino with MSSM fields (we also 
include the term f lT3|) worked out in the main text). Besides the relations already written in 
Section [2l we also use the Higgs decomposition 



= H+ cos /3, = H- sin (3, 
Hi = vi + ^ {H^ cos a - H^sma + iHl sin /3) , 

Hi = V2 + ^ {Hi sin a + H^ cos a + iHl cos 13) , (42) 
v2 

the sfermion rotation between the mass eigensates /i, /2 and the interaction eigenstates 
and the relations between the gaugino, and neutralino/chargino mass eigenstates [37j 



B = {N*,PL + KaPn)Xi 
= {N:,Pl + N,2Pr)x- 



A^ + lA^ 

We find 



= w={v;,p, + u,,Pn)x, 

= W'={u;,Pl + V,iPr)x- (44) 



Ant = ^Wx~^ ^jx" ~^ ^Zx° ~^ ^Hx~^ ^HOx° ~^ ^ff ~^ ^GG (^^) 
+^WHOx + ^WHx° + ^WWx° + ^Wff + -^T-H^X + ^-yff ~^ ^ZHx + ^ZH^''x° ~^ ^Zff 
+Cwyx + -^VFZx + ^GGG + '^G// 



^We denote the superfields associated to SM l.h. fermions by (PlJ^, t'L), {Phe, cl), (Plu, ul), {Phd^ dj^ , 
and the superfields associated to SM r.h. fermions by (Pte^, ejj), {Plu"^, u^), {Phd'^, d*j^; family indices 



are understood. 
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where 



' Mp 



(V2Mw cos ^ U;,r)^^ + J [Y: r]ru;,dp] Pj 



+ [y2Mw sin/3 V.^^"^ + \ [7", rh'^V^id^^^^ Pr 
+ f v^M^ sin /3 V;, + ^ [7^ 7'^] r V*,dp) pJ 1^" Xj 



M, 



(46) 



^Tx° = -TJT- t^"' + ^i's) + (cos 9^ Na + sin 9^ iV^a) Pi?] x° 

4 ikZ p 



(47) 



• -^ij, { [Mz (cos PN*,- sin ^ TV* ) rj'^'^ 



+7 [7", 7l 7'' (cos 9^ N*^ - sin 9^ TV* ) af )] P^ + [Mz (cos /3 TV.3 - sin P N,^) rj^'^ 



+1 [Y, r]Y {cos9^N,2-sm9^Na)df^]PR]z„ 



~(o) 



(48) 



'^^^ = -^i/;^v^i9'^//+(sin^C/;2Pi-cos^y,-2Pii)x,^ + h.c. 



(49) 



J^HOx° = - irr « ^i'a + sin « A^a) - (cos a iV,3 + sin « N^^) Pr] d^'Hl x ■ 

+i [(- sin a A^*3 + cos a N*^) Pj, - ( - sin a + cos a A^j4) Pr] ^''if^ Xi 
+ [(sin ^ TV* + cos ^ TV*4) Pl + (sin ^ TV.3 + cos ^ TV,4) Pr] d'^Hl x° } (50) 



5^ {cfd^f* + Sfd'^f;)p^f+ J2 {-Sfd'h + Cfd^h)PLr 



J=v,e,u,d 

c 



f=e,u,t 



+h.c. 



GG 



A Ml 



(51) 
(52) 



C 



WH°x 



• ^ (cos a U* Pl + sin a Vj2 Pr) + {- sin a U* Pl + cos a Vj2 Pr) 

+i (- sin ^ [7*2 + cos P Vj2 Pr) H^] x ■ + h. c. (53) 
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1 



WHx^ - 

£,wwx° — 



Mp 



g (cos 13N*^Pl + sin /3 iV,3 Pr) W''^ H' x° + h. c. 
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■ [7^ 7l 7'' (A^i^ + Pn) W; W; x° 



A Ml 



(54) 
(55) 



-i- ^ I [(c. i>i* + i>2*) Pl e + (c„ {it + s^u*^) Pl d\ 



(ce e\ + Se e^) Pl + (ca dl + Pl } + h- c. (56) 



■'iHx 



Mf 



ij^ V2 e (- sin /3 U*^ Pl + cos /3 Vj2 Pr) H+ Xj + h- c. 



(57) 



-iff 



-ZHx 



-^^|},V2eA^ J2 Qf [{cf fl + Sf f;) Pl f - {sf f\ + Cf h) Pl f] +h.c. 



f=e,u,d 

1 - g (cos^ 9yj — sin^ 9yj) 



(58) 

(- sin /3 [7*2 Pi, + cos /3 l^.a P^) H+ x^j + h. c. (59) 



CzH°x° 



Mf 



V2 cos 9yj 

o ^ I [(cos a 7V;3 - sin a N*^) Pl + (cos a N^s - sin a N,^) Pr] x • 

Z cos C/^y 1. 

- [(sin a N*^ + cos a N*^) Pl + (sin a Nis + cos a Ni^) Pr] Z*" x° 
-I [(sin /3 iV^, - cos /3 TVJ^) P^ - (sin /3 7V,3 - cos /3 TV^) Pr] ^'^ ^^3° X ° } 

(60) 



'Zff 



1 - \f2gZ^' 



^WZx — 



COS ^„ 



E zi{cff:+sff;)PLf- zi{-sfh+cff,)PLr 



.f=v,e,u,d 



f=e,u,d 



+h.c. 
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[7", 7l 1^ sin ([//i Pl + l^i + h-c 



4Mp 
9 



GGG 



29s_ 

8 Ml 



[Y, 7l7T''G«GjG^ 



(61) 
(62) 

(63) 
(64) 



Gff 



Mp 



^G-^ J2 [ {cf A"* A* + Sf A- /*) ^PLf--^-SfX- h + CfXy~2), Pl f^] +h.c. 



f=u,d 



(65) 



In the above expressions, d^^\ and d^^^ denote a derivative acting only on the W, and Z 

/ _ 1 _ 

Isin^^,, -i ' i^;v.2fl f_ . ^ 7/ 



fields, respectively; qj denotes the electric charge of the fermion /; Zj^ = |, — | + sin^ duj, | 

„ -fsin^^,,,. isin^^^ for 



y-w, 2 ' 3 



sin 9u, for v. e, u. d, respectively; Z 



R 



sm 



e, M, d, respectively; A" are the Gell-Mann matrices. Moreover, in some of the above ex- 
pressions we have also used the identities 

(^M Pl Xjy = Xj Pr = Pr X • , (^m [7", 7^] Y Pl X •) ^ = -^^ b", 7^] Y Pr Xj 

(66) 
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B Explicit derivation of the effective goldstino-matter 
Lagrangian in the non- derivative form 



Statement. 



9£» 



d, 



dC 



MSSM 



— does not contribute to S-matrix elements, at all orders 



in perturbation theory (with the only restriction that no goldstino enters in propagators), for 
arbitrary initial and final state, with one goldstino external line. 



Proof. Specifically, we need to show that 



(/|T exp 



d xC 



int. 



dC 



MSSM 



I I A 

— xMv 



(67) 



where £mssm = -^free + -^int., e is the global SUSY variation parameter, $j denotes any of the 
MSSM fields, and x is the goldstino. This proof is necessary to go from eq. to eq. fl25l) 
in the main text. The term with £free includes the free action for the MSSM fields. 
We have 



dC 



MSSM 



dC 



MSSM 



9^ free rj dUfj-gg 



free e.o.m. of $i + 



+ 

dCint. 



dC 



int. 



9$, 



dC 



int. 



9„ 



dC 



int. 



(68) 



The operator produced by the variation of £free is the (classical) equation of motion for the 
free field For this reason, we denoted it as "free e.o.m. of We actually prove that 



(/IT 



{,[./ 

-(/iTlexp 



d xC 



int. 



d^y [ free e.o.m. of $j] ^—^x \ 



dS 



int. 



9$, 



dC 



int. 



—x}\i), 



(69) 



from which eq. fl67|) immediately follows. We work out the l.h.s. of this expression, and 
we show that it is equal to the r.h.s.. When we use Wick's theorem to eliminate the time 
order product, the operator "free e.o.m. of $j" either acts on the initial or final state, or 
it is contracted with the field $j inside exp J d^xCint.], present in the exponent. In the 
former case, one obtains zero, since the fields in the initial and final states are free fields, 
whose wave functions obey the free equations of motion. The contraction gives instead a 
nonvanishing contribution. We note that "free e.o.m. of is linear in <I>j, so only one 
contraction with a single term in exp [i J d^xCi^t] takes place. We therefore have (normal 
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ordering is understood) 



l.h.s.of® ^ (/|T^xp[./.^.AL.]/A[freee.o.n.. oU.] M^) 



(/|t|5^^[^ y d'xC^^Y j A [free e.o.m. of $,] 
(/iTjf^f^i^^^^^^nz / d'zcL I d'y [ free e.o.m. of —X ) 



. n=l 



(70) 



where in the last step we have used the fact that "free e.o.m. of $i" contracts with all the n 
actions appearing in the n— th term in the expansion series of the exponent. Since the sum 
in the last expression is again the expansion series of the exponent, we have found that 



l.h.s. of dSnD = exp 



i J d'^xCint. i J d'^zCiat. J d'^y[ free e.o.m. of $i] ^—^x \ K) 



(71) 

To proceed, we need to recall the specific dependence of the interaction Lagrangian on 



For MSSM fields, we have 



(72) 



where n is an integer, and where the coefficients A and can depend on fields other than 
$i. Since the operator "free e.o.m. of $j" does not contain any other field rather than $j, 
these coefficients do not participate to the contraction. Notice that in the last step we have 
disregarded a boundary term that does not contribute to the interaction action. Under the 
time ordering, we therefore have 



J d'^zCint. (z) J d^y [ free e.o.m. of $j 



d^z 



Y,An {z)n<!>r' {z)-d^B^{z) 



{z) j d^y [ free e.o.m. of (y)] (73) 



We can now use the fact that [$i(^) free e.o.m. of ^i{y)] = i5^^^ (y — z). This is imme- 
diate from the fact that the contraction of ^i{z) with the operator ^i{y) entering in the 
expression in square parenthesis is the propagator, which is the "inverse" of the operator 
that forms the equation of motion (for instance, □ + m^, if $j is a scalar). We also use the 
fact that 



Y,An{z)n^r'{^)-d,B^{z) 



int. 



9(9^$. 



(74) 
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as can be immediately seen from eq. ( I72|) . Therefore, eq. ( !73|l can be continued to give 
I } I 



j d'^zCint. (z) J d^y [ free e.o.m. of $j 
Inserting this into eq. fl7T]) . we finally have 



l.h.s. of dSni) = {f\T |exp 

= r.h.s. of 
This completes the proof. 



i / d'^xC 



int. 



d yi 



- d„ 



(75) 



i / d'^yi 



o 



— X r 



(76) 



C Amplitudes for W — )► gravitino + neutralino 

The amplitude for the scattering {k) + W~ {k') G {p) + x° (g) can be written in the 
form 

iMi = % [Mi^i-t + M^,lz + Mi,2 + -M*,3 + -Mm + M^,^H^ + A1i,5//2] (77) 

The index i on the total and the partial amplitudes denotes the neutralino mass eigenstate 
produced in the reaction. The numerical index on the partial amplitudes denotes the order 
of the corresponding diagram in Figure 121 Notice that the first and last diagram in the 
Figure correspond to two different terms in ( 1771) . 

In writing the partial amplitudes, we treat Majorana spinors as explained in [IQ]. The 
Feynman rules for the vertices with the gravitino are immediately obtained from the terms 
listed in Appendix [Aj after eq. ( H6|) . The MSSM vertices entering in the diagrams are instead 
obtained from |3] and P7] 



MSSM 



-ig cos Q^Z- [2d^Wt - d,W^ - n^^d^W^) - W^^ (29^1^- - d^W;^ - v^.d^W 
-9 W- [Og Pl + Pr] rx'^ + g W; % r [Of; Pl + Of* Pr] x° 
+g Mw W+ W^- [cos (/3 - a) + sin (/3 - a) if °] 



where 



of^^-^N,,v;, + N,,v;, 



os^^n:,u,2+n:,u,. 



(79) 



and we remind the reader that the matrices A^, U, V are defined in eqs. 
We find 



iM 



4Mf 



[(cos N*^ + sin N*^) Pl + (cos 6^ + sin 9^ A^^a) Pb] Vi (g) 



{p + qY 



-% {2k + k')^ bi + « (A; + 2A;')„ < + ^ (A; - k!)" r]^p (k) {k') (80) 
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iz 



Mz (cos - sin fiN*,) r]^'^ + + Y] Y (cos O^N*, - sin O^N*,] 



+ 

-9aa' + 



Mz (cos pNis - sin pN^) + | + ^, t''] t'' (cos Q^N.^ - sin A^, 



Pl 
Pr 



■g cos \-i {2k + k')^ (5^' + i{k + 2k') ^ 5"^' +i{k- k'f rj^p] (k) {k') (81) 



iMi 



i,2 



Mr 



7a] 7^ 



(^72 Mm. cos /3U*^ 6^ + ^[f(:, 
+ (^^/2 M^^ sin (3 Vj2 S^ + ^[$, 7a] 7^^ ^ii) ^i? 



(P - kf - Ml 



-ig) [0§ Pl + 05 Pk] 7^ v, (q) e% (k) {k') (82) 



iMi 



i,3 



Mf 



V2 Mh. cos /3 Uj2 + - [H\ 7^] 7^^ Uji ) ^ii 



+ Mw sin ^ y.* 5^ + 1 [H\ 7/3] 7^^ 1^/1) 
— U' + M- ^ 

(^^) 7" [Og* Pl + Oj* Pn] Vi (q) (k) et {k') (83) 



{p-k'Y -M\ 



Xj 



^•^M = V'/.(P) [7a,7^]7''(iV;^i^L + iV,2P«)^i(g)e°(A;)e^(A;') (84) 



4M, 



« A4i,5ifi = {p) TT- [(cos a iV;3 + sin a N*^) Pl - (cos a + sin a Ni4^) Pr] {-i) {k + k'/ 
IVlp 



■Vi (q) 



{k + k'Y - Ml, 



i 9 Mw gap cos {p - a) e\ (k) et {k') 



(85) 



i Mi,5H2 = 'iPn (p) -ry- [(- sin a N*^ + cos a N*^) Pl - (- sin a Ni^ + cos a Ni^) Pr] (-i) {k + k'f 
IVlp 



■Vi (q) 



{k + k'Y - Ml, 



i g Mw gap sin (/3 - a) e° {k) i_ {k') 



(86) 



where t/'^, Vi, e", e^! denote the gravitino vector spinor, the neutrahno spinors, and the 
and W~ polarization vectors, respectively (the spin and polarization indices are understood). 
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